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ABSTRACT 

We present an empirical s-process abundance distribution derived with explicit knowledge of the 
r- process component in the low-metallicity globular cluster M22. We have obtained high- resolution, 
high signal-to-noise spectra for 6 red giants in M22 using the MIKE spectrograph on the Magellan- 
Clay Telescope at Las Campanas Observatory. In each star we derive abundances for 44 species of 
40 elements, including 24 elements heavier than zinc (Z — 30) produced by neutron-capture reactions. 
Previous studies determined that 3 of these stars (the "r + s group") have an enhancement of s- 
process material relative to the other 3 stars (the "r-only group" ) . We confirm that the r + s group 
is moderately enriched in Pb relative to the r-only group. Both groups of stars were born with the 
same amount of r-process material, but s-process material was also present in the gas from which the 
r + s group formed. The s-process abundances are inconsistent with predictions for AGB stars with 
M < 3 M Q and suggest an origin in more massive AGB stars capable of activating the 22 Ne(a,rt) 25 Mg 
reaction. We calculate the s-process "residual" by subtracting the r-process pattern in the r-only group 
from the abundances in the r + s group. In contrast to previous r- and s-process decompositions, 
this approach makes no assumptions about the r- and s-process distributions in the solar system and 
provides a unique opportunity to explore s-process yields in a metal-poor environment. 
Subject headings: globular clusters: individual (NGC 6656) — nuclear reactions, nucleosynthesis, 
abundances — stars: abundances — stars: AGB and post- AGB — stars: Popula- 
tion II 



1. INTRODUCTION 

By the time the average global star formation rate 
peaked in galaxies destined to grow to the size of the 
Milky Way 1-3 Gyr after the Big Bang, vigorous heavy 
metal enrichment had already begun. Elements heav- 
ier than the Fe-group are traditionally understood to be 
produced mainly by two processes, the rapid and slow 
neutron-capture processes. The r-process ("r" for rapid) 
manufactures heavy nuclei by overwhelming existing nu- 
clei with a rapid neutron burst on timescales ~ 1 s, far 
shorter than the average /3-decay timescales that could 
return unstable nuclei to stable ones. The s-process ("s" 
for slow) manufactures heavy nuclei by adding neutrons 
to existing nuclei on timescales slow relative to the aver- 
age /3-decay rates. Each of these two neutron (n) capture 
processes contributes about half of the heavy elements 
in the solar system (S.S.), which samples the chemistry 
of the interstellar medium (ISM) at one point in the 
Milky Way disc more than 9 Gyr after the big bang. 
The r-process requires an explosive, neutron-rich envi- 
ronment, suggesting an association with the core col- 
lapse supernovae (SNe) that claim the lives of massive 
stars (M > 8 M ), while the s-process may be acti- 
vated in less massive stars (1 < M < 8 Mq) during their 
asymptotic giant branch (AGB) phase of evolution. En- 
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richmcnt of the ISM by r-process material may begin a 
few tens of Myr after star formation commences, while 
s-process enrichment requires at least 50 Myr to several 
Gyr depending on the AGB mass ranges involved . 

Afte r an early description of the s-process by [ B 2 FH 
(fil^ . lClavton et all ((1961 and lSeeeer et al.1 (|1965D de- 
veloped the phenomenological (also known as the "clas- 
sical") approach that dominated s-process modeling for 
decades to follow. This method takes advantage of the 
fact that the product of the n-capture cross section and 
the s-process abundance of each isotope is slowly variable 
and can be approximated locally as a constant. These 
authors also recognized that a single neutron flux is in- 
suffi cient to reproduce the s-onl y isotopes in the S.S. (See 
also lClavton fc Rassbachl[i967l ) In order to explain the 
s-process distribution in the S.S., at least 3 components 
are required, known today as the "main," "weak," and 
"strong" components. The main component accounts 
for isotopes from 90 < A < 207, the weak component 
accounts for the bulk of the production of isotopes with 
A < 90, and the strong component accounts for more 
than half of 208 Pb. 

More and improved experimental data collected 
over subsequent decades revealed the shortcom- 
ings of this phenomenologic al approach (e.g., 
Kap peler. Beer, fc Wisshaki 119891 ). Predictions of 
s-process yields made by post-processi ng stellar evolu- 
tion m odels with reaction networks (e.g.. lArlandini et al.l 
1999) improved the fit, particularly near the closed neu- 
tron shells at N = 50, 82, and 126. Eventually the 
full reaction networks w ere integrated into the stellar 
evolut ion codes (e.g., iStraniero. Gallino. fc Cristallol 
120061 : iCristallo et al.l 120091 ). Nucleosynthesis via the 
s-process depends on a number of variables including 
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mass, metallicity, and s-process efficiency. Uncertainties 
in the mass dredged up after each thermal instability 
(which brings s-process material to the surface) and the 
mass-loss rate further complicate predictions. Detailed 
models are constrained by spectroscopic observations of 
s-process material in intrinsic (i.e., self-enriched) stars 
or extrinsic (i.e., enriched by a binary companion or 
born with the s-process material) ones. These models 
have mainly focused on low- and intermediate-mass 
AGB stars (i.e., < 3 M Q ) and are quite successful at 
reproducing both the S.S. s-process isotopic distribution 
and the el emental distributions observed in a variety of 
stars (e.g.jSmith fc Lambert! 1 19861; iLambert et alj|1995t 
IBusso et alJ Tl999. 2001; Bist erzo et al. 20091|2011|). 

The s-process efficiency is largely governed by the 
conditions that activate reactions to liberate neutrons, 
and two sources of neutrons have been identified in 
AGB stars. The first, the 13 C(a,n) 16 reaction, is ac- 
tivated at temperatures around 1x10 s K. 13 C is of pri- 
mary origin, synthesized from proton captures on freshly 
produced 12 C. The 13 C pocket is thought to form in 
the top layers of the region between the H and He 
shell-burning regions when protons from the H enve- 
lope are mixed into this region during the third dredge 
up. The amount of 13 C in the pocket can be thought 
of as one measure of the s-process efficiency. The 
other source, the 22 Ne(a,n) 25 Mg reaction, is activated 
at somewhat higher temperatures near 3.5 xlO 8 K. 22 Ne 
is also primary. It is produced by the reaction sequence 
14 N(a,7) 18 F(£+zy) 18 0(a,7) 22 Ne, where 14 N is also pri- 
mary as the most abundant product of CNO burning. 
Neutron densities from the 13 C and 22 Ne reactions may 
reach ~ 10 7 cm" 3 and > 10 11 cm" 3 , respective ly . See , 
e.g., reviews by IBusso. Gallino. fe W asscrburgl (|1999f ) 
and iStraniero et aTT l 20061 ) for further details. 

The heavy elements in the S.S. are the products of 
many and various stars, and the stellar evolution param- 
eter space necessary to fully reproduce the S.S. s-process 
pattern is vast and gradually being explored. Only in 
the S.S. is the complete heavy element inventor y known 
with great precision at the isotopic level (e.g.. lLoddersI 
2003). Isotopes that can only be formed by the r- or 
s-process are readily identified, but no element in the 
S.S. with 30 < Z < 83 owes its presence entirely to 
the r- or s-process. Limited by the assumption that 
only two processes contribute, the r- and s-process con- 
tent in S.S. material can be estimated by the formula 
N Qy r = iV Qi tot — Nq iS . That is, the r-process "residual" 
equals the total S.S. abundance minus the s-process con- 
tribution, which is ob tained by either phenomenological 
or stellar models (e.g . , iSeeger et al.lll965l lCameror]|1973l 
iKappeler et al.lll989l ). Nearly all abundance information 
in other stars is in the form of elemental abundances. In 
certain astrophysical environments only one process or 
the other contributes, enabling direct comparison with 
model predictions. The difficulty lies in identifying suit- 
able stars whose heavy elements may be reliably inter- 
preted as having originated in only one process or the 
other. 

One such star, CS 22892-052, with a metallicity less 
than 1/1000 solar ([Fe/H] = -3.1)B was di sc overe d 
in the survey of iBeers. Preston, fc Shectmanl (|1992j) . 

5 We adopt standard definitions of elemental abundances and 



CS 22892-052 has a heavy element abundance pattern 
that very nearly matches the scaled r-process residual s 
in the S.S. (e.g.. ISneden"et~aT1[l99l [Cowan et allll99ll) . 
Several other metal-poor stars with this pattern have 
been found, and nearly all stars analyzed to date con- 
tain detectable quantities of elements heavier than the 
Fe-group. These elements are frequently attributed to 
r-pro c ess nucleosynthesis fe.g.. | Trura n 1981, Mc William 
1 19981 ISneden, Cowan, fc Gallino! 120081 IRoederer et al. 
2010b). The consistent r-process abundance pattern ob- 
served in several stars heavily enriched by r-process ma- 
terial inspired the idea that r-process abundances every- 
where (at least for Z > 56) may be scaled versions of 
the same pattern; however, stars with less extreme levels 
of r- process enrichment clearly deviate from th is pattern 
(e.g., iHonda et al1l2007l IRoederer et~aDl2010bf ). 

Some metal-poor stars contain s-process material 
mixed with the r-process contribution. Obtaining an em- 
pirical measure of the s-process content in stars other 
than the sun is difficult because a level of r-process 
enrichment must be assumed. The metal-poor glob- 
ular cluster (GC) M22 provides an opportunity to 
probe s-process enrichment in a low-metallicity environ- 
ment where the r-process content is explicitly known. 
Recent spectroscopic studies have demonstrated that 
star-to-star variations in h eavy elements exist in M22 
(|Marino et all 120091 l2011bl) . This metal-poor ([Fe/H] = 
— 1.76 ± 0.10) GC hosts two groups of stars, each with 
different amounts of heavy elements (Y, Zr, Ba, La, Nd) 
that i n the S.S. are overwh el mingly due to the s-process 
(e.g., iSimmerer et al.l [2004T ). iMarino et al.l showed that 
the abundances of these elements, together with the to- 
tal C+N+O and overall Fe-group abundances, increase 
as a function of metallicity. In contrast, [Eu/Fe] has no 
metallicity dependence (only 3% of S.S. Eu was produced 
by the s-process) , demonstrating that the heavy element 
variations are due to different amounts of s-process ma- 
terial. Thus, the chemistry of M22 suggests one stellar 
group formed from gas enriched by r-process nucleosyn- 
thesis and a second group formed from gas also enriched 
in s-process material. Multiple stellar groups in M22 are 
also revealed in a split in the sub-giant branch (SGB) 
revealed by Hubble Spac e Telescope photometry (|Piottol 
2009: IMarino et al.ll2009h . 

The chemical pattern revealed in M22 makes this GC 
a suitable target to investigate s-process abundance dis- 
tributions. Observations indicate that the r-process con- 
tent of both stellar groups in M22 is the same. Observa- 
tions also indicate that the more metal-rich group (here- 
after referred to as the "r + s group" ) was formed from 
gas also enriched by s-process material. We can subtract 
the r-process abundance pattern (established empirically 
in the metal-poor group, hereafter referred to as the "r- 
only group") from the abundance pattern in the r + s 
group to derive an empirical s-process abundance distri- 
bution. One favorable aspect of this approach is that it 
does not rely on the decomposition of S.S. material into r- 
oi s-process fractions to interpret abundances elsewhere 
in the Galaxy. 

2. OBSERVATIONS 



ratios. For element X, loge(X) = log 10 (Nx/N u )+ 12.0. For ele- 
ments X and Y, [X/Y] = log 10 (7V x /iVY)* - log lo (JV x /iV Y ) . 
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TABLE 1 

Photometry, Atmospheric Parameters, Radial Velocities, Exposure Times, and S/N Estimates 



Star 


V 


(B - V) 


T e g 


log 9 


«t 


[Fe/H] 


RV 


^cxp 


S/N 


S/N 


S/N 


S/N 








(K) 




(km s ) 




(km s _1 ) 


(s) 


(3950A) 


(4550A) 


(5200A) 


(6750A) 


1-27 


12.39 


1.28 


4455 


1.45 


1.60 


-1.73 


-127.8 


2600 


40/1 


95/1 


95/1 


210/1 


1-37 


12.01 


1.45 


4370 


1.05 


1.50 


-1.73 


-157.7 


3800 


50/1 


125/1 


140/1 


340/1 


1-53 


12.69 


1.36 


4500 


1.35 


1.55 


-1.74 


-145.4 


3000 


45/1 


105/1 


120/1 


270/1 


1-80 


12.53 


1.38 


4460 


1.15 


1.55 


-1.70 


-149.8 


3100 


40/1 


90/1 


100/1 


230/1 


111-33 


12.25 


1.40 


4430 


1.05 


1.70 


-1.78 


-145.8 


1600 


45/1 


105/1 


120/1 


275/1 


IV-59 


11.93 


1.45 


4400 


1.00 


1.70 


-1.77 


-152.8 


1600 


45/1 


110/1 


125/1 


300/1 



Six probable members of M22 were observed with the 
Magellan Inamori Kyo cera Echelle (MIKE) spectrograph 
(|Bernstein et al.H2003D on the 6.5 m Magellan- Clay Tele- 
scope at Las Campanas Observatory on 2011 March 17- 
18. These spectra were taken with the 0.7" x 5.0" slit 
yielding a resolving power of R ~ 41,000 in the blue and 
R ~ 35,000 in the red, split by a dichroic around 4950A. 
This setup provides complete wavelength coverage from 
3350-9150A, though in practice we only make use of the 
region from 3690 to 7800A where the lines of interest are 
located. Data reduction, extraction, and wavelength cal- 
ibration were performed using the MIKE data reductio n 
pipeline written by D. Kelson. (See also IKelsonl [2003) . 
Continuum normalization and order stitching were per- 
formed within the IRAF environment 

The six observed stars are all cool giants on the 
M22 red giant branch (RGB). Table [T] lists the pho- 
tometry from the Stetson database (priv. commu- 
nication, corrected for differential reddening as in 
iMarino et al.l I2011TJ) atmospheric parameters (adopted 
from lMarino et all see Section [3]), Heliocentric radial ve- 
locities (RV), exposure times, and signal-to- noise (S/N) 
estimates for our targets. We estimate the S/N based on 
Poisson statistics for the number of photons collected in 
the continuum. We measure the RV with respect to the 
ThAr lamp by cross correlating the echelle order contain- 
ing the Mg I b lines in each spectrum against a template 
using the fxcor task in IRAF. We create the template 
by measuring the wavelengths of unblended Fe I lines in 
this order in star IV-59, which has the highest S/N in a 
single exposure. We compute velocity corrections to the 
Heliocentric rest frame using the IRAF rvcorrect task. 
This method yie lds a total uncertainty of 0.8 km s _1 per 
observation (see iRoederer et al1l2010aD . Our RVs are in 
good ag reement (A = 0.7 ± 0.7 km s _1 ) with those de- 
rived by IMarino et al.1 (|2009f) for 4 stars in common. 

3. ABUNDANCE ANALYSIS 

We perform a standard abundance analysis on the 
6 stars observed with MIKE. We adopt th e atmo - 
spheric parameters derived by IMarino et al.1 (|2011b[ ) 
and use ct-enhanced AT LAS9 model atmospheres from 
ICastelli fe Kuruc3 (|2004l ). We perform the analysis us- 
ing the latest v ersion of the spectral analysis code MOOG 
(|Snedenl H973) . with updates to the calculation of the 
Rayleigh scattering contribution to the continuous opac- 

6 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



ity described in lSobeck et al.l (|2011| ). We measure equiv- 
alent widths (EWs) by fitting Voigt absorption line pro- 
files to the continuum-normalized spectra, and we derive 
abundances of Na I, Mg I, Al I, Si I, K I, Ca I, Ti I and 
II, Cr I and II, Fe I and II, Ni I, and Zn I from a stan- 
dard EW analysis. Abundances of all other elements are 
derived by spectral synthesis, comparing synthetic spec- 
tra to the observations. This is necessary for species 
whose lines may be blended, have broad hyperfme struc- 
ture (HFS), or have multiple isotopes whose electronic 
levels are shifted slightly. The linelist, atomic data and 
references, and derived abundances for each line are pre- 
sented in Table [21 which is available in the online-only 
edition. 

In Figure[Tjwe show synthetic spectra fits to the region 
around the Pb I line at 4057A. The solid lines represent 
the best-fit abundance, and the dashed lines represent 
variations in this fit by 0.3 dex. The 6 stars in our sam- 
ple have very similar atmospheric parameters. In our 
syntheses we adjust the linelist to fit blending features 
in one star and leave these adjustments unchanged in 
the analysis of other stars. This preserves a differential 
quality in the abundance analysis, which is important in 
the case of abundances derived from very few or heavily 
blended lines. 

The EWs measured by IMarino et al.1 (|2011b[ ) from a 
high S/N VLT/UVES spectrum of 1-27 are systemati- 
cally higher by 4.9 ± 0.6 mA (cr = 3.2 mA). the only 
case where the offset is larger than the standard devia- 
tion of the residuals (3.2-3.4 mA in all stars). In 1-27, 
this translates to no significant difference in the derived 
[X/Fe] ratios (since both abundances are similarly af- 
fected), A[X/Fe] = 0.00 ± 0.01 (a = 0.03 dex), and 
the change in metallicity is A [Fe/H] = 0.06. EWs mea- 
sured from the lo wer S/N APO/ ARC Echelle spectrum 
of 1-27 analyzed bv lMarino et al.l are systematically lower 
by 5.2 ± 2.6 mA (a = 11 mA). Differences at this level 
may be expected when analyzing spectra of various reso- 
lution and S/N, collected over many nights using instru- 
ments on different telescopes in different hemispheres, 
so we do not pursue the matter further. We also com- 
p are our derived abundance ratios with those presented 
in IMarino et al.l After accounting fo r the different sets o f 
log(g/) values, S.S. abundances ( see IMarino et aU feOOg). 
and line- by-line mean offsets (see Appendix), the abun- 
dance offsets for most [Fe/H] and [X/Fe] ratios can be 
immediately accounted for. Offsets in other species (Ti I, 
Cu I, Zn I, La II, and Nd n) are unexplained by these fac- 
tors but probably result from the S/N and small num bers 
of features available in the spectra of IMarino et al.1 For 
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Fig. 1. — Comparison of observed (open squares) and synthetic (lines) spectra around the Pb I 4057. 8A line. The left panels show the 3 
r-only stars, and the right panels show the 3 r + s stars. The solid colored line (left panels, blue; right panels, red) indicates the best-fit 
abundance, the dashed lines indicate variations in the best-fit abundance by 0.3 dex, and the solid black line indicates a synthesis with no 
Pb I present. 



TABLE 2 
Line-by-line Abundances 



Species 


A 


E.P. 


log(ff/) 


Ref. 


1-37 


111-33 


IV-59 


1-27 


1-53 


1-80 




(A) 


(eV) 


















Na I 


5682.63 


2.10 


-0.71 


1 


4.68 


4.09 


4.95 


4.69 


4.70 


5.10 


Na I 


5688.20 


2.10 


-0.45 


1 


4.81 


4.19 


5.05 


4.88 


4.81 


5.18 


Na I 


6154.23 


2.10 


-1.55 


1 


4.63 


4.07 


4.98 




4.63 


5.04 


Na I 


6160.75 


2.10 


-1.25 


1 


4.77 


4.13 


4.89 


4.68 


4.67 


5.13 



References. (1) Fuhr & Wicsc (200 9]); (2) IChang fe Taiigi (119901); (3) 

ILawler fe DakirJ 119891) . using HFS from iKurucz fc Bell (1995); (4) IBlackweil eFaLl 
Tl9*82*alIbT) . increased by 0.056 de x according toTG*revesslTeTal.l dl989!1: f5^IPiekermg^t**aLl 
020011) , wi th corrections g i ven inlPicker ine et a l.l (120021) ; ( 6) Whaling ct al. (1985), using 
HFS from Kurucz & Bell (1995); (7) Sobcck ct al. (2007); (8) Nilsson ct al. (2006); (9) 
Blackwcll- Whitehead & Bcrgcmann (2007), using HFS from Kurucz & Bell (1995); (10) 
booth et al.l (119841). using HFS from IKurucz fe Bell (Tt995T ) fill IQ'Brian et alj ll99H1: 
(12)[Mclcridcz & Barbuv ( 2009 ?) ; (13 ) iNitz et al.l l|1999l ), u sing HFS from IKur ucz & Bell 
(199 5]); (141 lUardon et al.l 119821). using H FS from IKurucz fe Bel l (1995); (15 ) 
Wickliffe & Lawlcr (1997a); (16) Biclski (1975), using HFS from Kurucz & Bell (1995); 
(17) Bicmont & Godcfroid (1980); (18) Migdalck & Baylis (1987); (19) Hannaford ct al. 
119821) ; (20) IBiemont et al.l 119811 1; (2 1) ILiung et al.| 12006?); (22 1 I Whaling fe Brault 
lT§§a) ; (23) [Wickliffe ct al. ( 1994:1; ( 241 IDuquette fe L awler (1985); (25) IFuhr fe Wiese 
120091). using HFS from IMcWilliarol 119981); (2 6) ILawler et al.l 12001al1 , using HFS from 
Ivans ct al. (2006); (27) Lawlcr ct al. (2009); (28) ILi et aTT fepOTI), using HFS from 
Sncdcn ct al. (2009); (29) Ivarsson ct al. ( 200IJ1, using HFS from Sned en et al] |[2009l; 
(30) IDen Hartog et all 12003 ); (31 ) ILawler eTa l. (2006); (32) Lawlc r et al.l 12001b|), u sing 
HFS from llvans et all 1200611; ( 33)IDen Hartog et alJ 120061); (34) Lawlcr et al. (2 0"01cD. 
using HFS from ILawler et al.l 12001dl) ; (3 5) I Wicklifte et alj 120001) ; (36) ILawler et al. 
1 2001) for bot h log(.g/1 value and H FS; (371 ILawler et al.l l|2008l) 7T38) Wickliffe fe Lawlei 
1 1997bl1 ; (3 91 ISneden et al.l (f2059h for both logfaf) value a nd H FS; (40) ILawler et ah 
l 20071) ; (411 llvarsson et al.l (120031) ; (421 IBiemont et al.l (f2000H ; (43l lNilsson et al.l 12002TT 

Note. — The full version of this table is available in the online edition of the journal. 
A limited version is given here to indicate the form and content of the data. Abundances 
are given as loge notation. A ":" indicates the derived abundance is less secure, and we 
estimate an uncertainty of 0.2 dex. See text for details. 
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M22 LINE -BY -LINE DIFFERENTIAL ANALYSIS 



TABLE 3 

Effect of Different Analysis Tools on Derived 
Abundances in 1-37 



A (A) 


loge A 


loge B 


log e c 


loge D 


La II Lines 


3988.51 


-0.74 


-0.66 


-0.71 


-0.61 


qooc: i a 
oyyo. 1 4 


— U. 1 U 


— U.DO 


— u.oy 


— U.ol 


4086.71 


-0.78 


-0.70 


-0.70 


-0.56 


4322.50 


-0.70 


-0.67 


-0.68 


-0.66 


4662.50 


-0.54 


-0.50 


-0.55 


-0.46 


4748.73 


-0.65 


-0.65 


-0.62 


-0.70 


4804.04 


-0.44 


-0.43 


-0.43 


-0.53 


4920.98 


-0.34 


-0.31 


-0.28 


-0.25 


4986.82 


-0.49 


-0.48 


-0.47 


-0.41 


5114.56 


-0.48 


-0.45 


-0.43 


-0.42 


5290.84 


-0.73 


-0.69 


-0.71 


-0.64 


5303.53 


-0.45 


-0.46 


-0.44 


-0.45 


6262.29 


-0.44 


-0.45 


-0.41 


-0.41 


6390.48 


-0.39 


-0.42 


-0.38 


-0.40 


Mean: 


-0.56 


-0.53 


-0.54 


-0.50 


cj: 


0.15 


0.12 


0.15 


0.12 


a/VN: 


0.040 


0.033 


0.038 


0.033 


Ce II Lines 


4073.47 


-0.35 


-0.24 


-0.32 


-0.23 


4083.22 


-0.25 


-0.17 


-0.21 


-0.14 


4120.83 


-0.20 


-0.17 


-0.17 


-0.08 


4127.36 


-0.24 


-0.20 


-0.24 


-0.15 


4137.65 


-0.26 


-0.18 


-0.20 


-0.13 


4222.60 


-0.23 


-0.19 


-0.20 


-0.16 


4364.65 


-0.21 


-0.18 


-0.19 


-0.15 


4418.78 


-0.23 


-0.15 


-0.16 


-0.08 


4486.91 


-0.26 


-0.24 


-0.24 


-0.23 


4560.96 


-0.21 


-0.19 


-0.18 


-0.17 


4562.36 


-0.17 


-0.14 


-0.14 


-0.08 


4572.28 


-0.16 


-0.11 


-0.02 


±0.00 


4582.50 


-0.05 


-0.02 


-0.05 


±0.03 


4628.16 


-0.08 


-0.07 


-0.02 


±0.04 


5274.23 


-0.20 


-0.21 


-0.22 


-0.20 


5330.56 


-0.22 


-0.24 


-0.19 


-0.22 


Mean: 


-0.21 


-0.17 


-0.17 


-0.12 


a: 


0.07 


0.06 


0.08 


0.09 


a/VN: 


0.017 


0.015 


0.020 


0.021 


loge(La/Cc): 


-0.35 


-0.36 


-0.36 


-0.38 


±0.043 


±0.036 


±0.044 


±0.039 


A MOOG with scattering 


, MARCS model 





B MOOG without scattering, MARCS model 
c MOOG with scattering, ATLAS9 model 
D MOOG without scattering, ATLAS9 model 



the purposes of the present study, we will focus on the 
internal abundance differences derived from our MIKE 
spectra. 

Table [3] shows the internal abundance precision pos- 
sible with this method when large numbers (N > 10) 
of lines are available across the visible spectral range. 
For this test we derive the abundances of La II and 
Ce II in star 1-37 using two grids of m odel atmo- 
spheres (MARCS IGustafsson et all l2008t ATLAS9, 
ICastelli fe KuTucl] 120041) and different treatments of 
Rayleigh scattering in MOOG. The elemental abun- 
dances and ratios are not dependent on the choice of 
model atmosphere grids or the treatment of Rayleigh 
scattering. For example, as shown in Tabled in response 
to the different analysis tools the derived log e(La/ Ce) ra- 
tio changes by no more than 0.03 dex, which is smaller 
than the statistical uncertainties (0.036 to 0.044 dex). 
Furthermore, the stars in our study were chosen to have 
similar colors (1.36 < (B — V)o < 1-52), metallicities 
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Fig. 2. — Differences in the mean abundances between the 3 
r-only stars and the 3 r ± s stars as a function of atomic number 
for K through Zn. Solid squares indicate neutral species, and open 
squares indicate singly-ionized species. The dotted line indicates 
zero difference. 

(—1.80 < [Fe/H] < —1.70), and atmospheric parameters 
(4370 < T eff < 4500 K and 1.00 < log g ^ 1.45; all based 
on the values presented in lMarino et al.|[2011bT )FI Thus 
a relative abundance analysis is appropriate, and in all 
subsequent discussion, tables, and figures we cite internal 
(i.e., observational) uncertainties only. 

Absolute uncertainties that account for errors in 
the derived a tmosph eric parameters are discussed in 
iMarino et~aLl (|2011b[ ) and presented in Table 4 of that 
work. In the present study, if only 1 line of a particular 
species has been measured we adopt an uncertainty of 
0.11 dex. This estimate is based on the mean standard 
deviation of individual lines for well-measured n-capture 
species (i.e., N > 3). Some lines in Table [5] are marked 
with ":" to indicate that the derived abundance is less 
certain due to significant blending features, difficult con- 
tinuum placement, etc. These lines have an adopted in- 
ternal uncertainty of 0.2 dex. 

One difficulty that is not minimized by our approach is 
that of comparing abundance ratios derived from species 
of different ionization states. Ratios of, e.g., [Ca/Fe] or 
[Eu/Fe] are computed by comparing Ca I to Fe I or Eu II 
to Fe II since both species of Fe are detected. Other ra- 
tios, such as [Pb/La] or [Pb/Eu], which compare Pb I to 
La II or Eu II, may be systematically uncertain. Note 
that for illustration purposes in the figures only we nor- 
malize the abundances of first-peak n-capturc elements 
observed in their neutral state (Sr I, Mo I, Ru I, and 
Rh i) to the singly-ionized abundances by the difference 
in Zr II and Zr I in each star (typically 0.2-0.4 dex). 

4. HEAVY ELEMENT ABUNDANCES IN M22 

In this section we analyze the abundance patterns in 
detail. The abundance results for each star in the r-only 
and r + s groups are given in Tables |4] and [5j respectively. 
Table[6]lists the mean abundances for each element in the 
r-only and r+s groups. We derive only upper limits from 
the Rb I line at 7800A. Due to blending by CN and CH, 
we are unable to derive abundances of Ir I or Th n in any 

7 For comparison, precision abundance analyses of nearby metal- 
rich dwarfs with stellar parameters similar to the sun often consider 
stars with T e g within 100 K, log g within 0.1 dex, an d [Fe/H] within 
0.1 d ex of the solar values to be "solar twins" (e.g.. iRamfrez et al.l 
2009). 
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TABLE 4 

Mean Abundances in the Three r-ONLY Stars 



Species 


Z 




1-37 








111-33 








IV-59 








(log e) 




1 V 


cr 


o> 






IV 


cr 


°> 


(log e) 


I V / L v. \ 


A 




cr 


o> 


Na I 


11 


4.72 


0.20 


4 


0.08 


0.041 


4.12 


-0.26 


4 


0.05 


0.026 


4.97 


0.53 


4 





.07 


0.033 


Mg i 


12 


6.25 


0.37 


3 


0.21 


0.121 


6.17 


0.43 


3 


0.17 


0.101 


6.16 


0.36 


3 





.20 


0.115 


A 1 T 

Al I 


13 


4.89 


0.16 


4 


0.09 


0.047 


4.50 


-0.09 


1 


0.11 


0.110 


5.33 


0.68 


1 





.11 


0.053 


hi I 


14 


5.89 


0.10 


3 


0.24 


0.140 


5.80 


0.16 


3 


0.17 


0.098 


5.86 


0.15 


3 





.18 


0.104 


K i 


19 


4.12 


0.81 


1 


0.11 


0.110 


3.93 


0.76 


1 


0.11 


0.110 


4.06 


0.83 


1 





.11 


0.110 


Ua I 


20 


5.04 


0.41 


8 


0.11 


0.038 


4.84 


0.36 


8 


0.15 


0.054 


4.93 


0.40 


8 





.10 


0.035 


C„ TT 

be II 


21 


1.68 


0.15 


5 


0.27 


0.120 


1.58 


0.20 


5 


0.30 


0.134 


1.59 


0.28 


5 





.29 


0.128 


Ti I 


22 


3.28 


0.05 


9 


0.07 


0.022 


3.11 


0.03 


9 


0.07 


0.022 


3.25 


0.11 


9 





.08 


0.027 


Ti II 


22 


3.79 


0.47 


9 


0.12 


0.039 


3.58 


0.40 


9 


0.10 


0.035 


3.64 


0.53 


9 


0. 


.09 


0.031 


V I 


23 


2.00 


-0.21 


5 


0.10 


0.045 


1.97 


-0.10 


5 


0.12 


0.052 


2.02 


-0.11 


5 





.09 


0.040 


Cr I 


24 


3.74 


-0.18 


6 


0.12 


0.049 


3.62 


-0.15 


6 


0.08 


0.033 


3.72 


-0.12 


6 





.07 


0.029 


Cr II 


2-1 


4.14 


0.13 


1 


0.11 


0.110 


4.00 


0.13 


1 


0.11 


0.110 


3.95 


0.15 


1 





.11 


0.110 


Mn I 


25 


3.13 


-0.59 


1 


0.12 


0.061 


3.12 


-0.45 


1 


0.12 


0.058 


3.16 


-0.47 


1 





.12 


0.061 


Fe I 


26 


5.78 


-1.72 


69 


0.13 


0.016 


5.63 


-1.87 


69 


0.10 


0.012 


5.70 


-1.80 


69 





.12 


0.015 


Fe II 


26 


5.87 


-1.63 


7 


0.17 


0.062 


5.73 


-1.77 


7 


0.08 


0.031 


5.66 


-1.84 


7 





.09 


0.034 


Co I 


27 


3.11 


-0.16 


4 


0.09 


0.043 


3.01 


-0.12 


1 


0.13 


0.065 


3.06 


-0.13 


1 





.18 


0.089 


IN 1 I 


28 


4.28 


-0.23 


10 


0.10 


0.033 


4.20 


-0.16 


10 


0.12 


0.039 


4.22 


-0.19 


10 





.11 


0.034 


tu I 


29 


1.83 


-0.65 


2 


0.08 


0.057 


1.64 


-0.69 


2 


0.08 


0.057 


1.74 


-0.65 


2 





.08 


0.057 


Zn I 


30 


2.89 


0.04 


2 


0.08 


0.055 


2.75 


0.05 


2 


0.08 


0.057 


2.80 


0.04 


2 





.08 


0.057 


Kb I 


37 


<1.00 


<0.20 


1 






<1.10 


<0.44 


1 






<1.10 


<0.38 


1 








Sr I 


38 


0.60 


-0.55 


1 


0.11 


0.110 


0.51 


-0.50 


1 


0.11 


0.110 


0.54 


-0.53 


1 





.11 


0.110 


Y II 


39 


0.34 


-0.25 


7 


0.04 


0.016 


0.18 


-0.26 


7 


0.07 


0.026 


0.18 


-0.19 


7 





.08 


0.030 


z^r I 


40 


0.69 


-0.18 


2 


0.08 


0.057 


0.78 


0.06 


2 


0.08 


0.057 


0.63 


-0.15 


2 





.12 


0.085 


Zr II 


40 


1.06 


0.11 


3 


0.09 


0.054 


1.04 


0.23 


3 


0.09 


0.050 


1.05 


0.31 


3 





.09 


0.050 


Mo I 


42 


0.33 


0.16 


2 


0.14 


0.099 


0.00 


-0.02 


1 


0.20 


0.200 


0.12 


0.04 


3 





.28 


0.164 


Ru I 


41 


0.09 


0.06 


3 


0.14 


0.082 


0.07 


0.18 


3 


0.13 


0.075 


0.18 


0.23 


1 


0. 


.20 


0.200 


DU t 
1X11 I 


45 


—0.84 


-0.18 


1 


0.11 


0.110 













—1.16 


-0.42 


1 


0. 


. 1 1 


0.110 


Ba II 


56 


0.76 


0.20 


2 


0.08 


0.057 


0.53 


0.12 


2 


0.15 


0.105 


0.53 


0.19 


2 





.18 


0.130 


La II 


57 


—0.54 


-0.01 


14 


0.14 


0.038 


—0.59 


0.09 


11 


0.08 


0.022 


—0.61 


0.13 


14 


0. 


.11 


0.029 


II 


58 


—0.17 


-0.13 


16 


0.08 


0.020 


—0.26 


-0.06 


16 


0.05 


0.013 


—0.29 


-0.03 


16 





.0 ( 


0.018 


Pr II 


59 


—0.87 


0.04 


4 


0.05 


0.025 


—0.87 


0.18 


4 


0.08 


0.040 


—0.93 


0.19 


1 





.05 


0.025 


Nd II 


60 


—0.13 


0.07 


24 


0.08 


0.017 


—0.18 


0.18 


21 


0.07 


0.014 


—0.22 


0.20 


24 





.08 


0.016 


O™ TT 

Ml II 


62 


—0.48 


0.19 


9 


0.07 


0.022 


—0.51 


0.30 


9 


0.07 


0.022 


—0.55 


0.33 


9 





.0 ( 


0.023 


Eu II 


63 


—0.85 


0.26 


3 


0.20 


0.113 


—0.88 


0.37 


3 


0.15 


0.085 


—0.92 


0.40 


3 





.17 


0.096 


Gd ii 


64 


-0.30 


0.25 


3 


0.14 


0.081 


-0.35 


0.35 


3 


0.16 


0.092 


-0.41 


0.36 


3 





.11 


0.061 


Tb ii 


65 


-1.48 


-0.15 


1 


0.11 


0.110 


-1.35 


0.12 


1 


0.11 


0.110 


-1.67 


-0.13 


1 





.11 


0.110 


Dy ii 


66 


-0.20 


0.33 


4 


0.13 


0.063 


-0.19 


0.48 


4 


0.09 


0.046 


-0.37 


0.37 


1 





.29 


0.144 


Ho ii 


67 


-1.10 


0.05 


1 


0.20 


0.200 


-1.15 


0.14 


1 


0.20 


0.200 


-1.20 


0.16 


1 





.20 


0.200 


Er ii 


68 


-0.30 


0.41 


2 


0.08 


0.057 


-0.48 


0.37 


2 


0.10 


0.070 


-0.55 


0.37 


2 





.07 


0.050 


Tm ii 


69 


-1.43 


0.10 


2 


0.18 


0.125 


-1.53 


0.15 


2 


0.25 


0.175 


-1.65 


0.09 


2 





.07 


0.050 


Yb ii 


70 


-1.05 


-0.34 


1 


0.20 


0.200 


-0.75 


0.10 


1 


0.20 


0.200 


-1.00 


-0.08 


1 





.20 


0.200 


Hf ii 


72 


-1.00 


-0.22 


1 


0.11 


0.110 


-0.93 


-0.01 


1 


0.11 


0.110 


-1.09 


-0.10 


1 





.11 


0.110 


Ir I 


77 


0.00 


0.34 


1 


0.11 


0.110 


-0.05 


0.43 


1 


0.11 


0.110 


0.15 


0.57 


1 





.11 


0.110 


Pb I 


82 


0.05 


-0.27 


1 


0.11 


0.110 


-0.01 


-0.19 


1 


0.11 


0.110 


-0.08 


-0.32 


1 





.11 


0.110 


Th ii 


90 


-1.55 


0.02 


1 


0.11 


0.110 


-1.42 


0.29 


1 


0.11 


0.110 


-1.46 


0.32 


1 





.11 


0.110 



Note. — Quoted uncertainties represent internal uncertainties only. ([Fe/H]} is listed in the ([X/Fe]} column for Fe I and Fe II. 

analysis. The differential results are listed in Table [7] and 
illustrated in Figure [5J When considering the standard 
error (a^ = a/^/N) of the mean line-by-line differen- 
tial abundances (column 5), [K/Fe], [Ca/Fe], [Sc/Fe]0 
and [V/Fe] show slight but significant (0.06-0.10 dex) 
enhancements in the r + s group, [Ti/Fe] and [Cr/Fe] are 
indistinguishable in the two groups, and [Mn/Fe] shows 
a slight (0.06 dex) defic iency in the r + s group. 

IMarino et~aT1 (|2011bH saw an increase by 0.06-0.15 dex 
in the [Cu/Fe] and [Zn/Fe] ratios of the r+s group, which 
we also detect. Furthermore, we find a similar — though 
smaller — enhancement in [Co/Fe] and [Ni/Fe]. The re- 
sults from a line-by-line differential analysis of these ele- 
ments are also listed in Table [7] and shown in Figure HJ 
A slight abundance enhancement in the Fe-group el- 

8 The Sc II lines give discordant abundances, which may indicate 
relatively large uncertainties in the log(gf) values, but the line-by- 
line results are extremely consistent. 



star in the r + s group. 

4.1. The Light and Fe-group Abundance Patterns in 

M22 

We derive abundances of Na I, Mg I, Al I, and Si I 
in each star of the sample. After accounting for dif- 
ferences in the log(q f) values between this study and 
IMarino et al.l (|2011b| ). these ab undances are in agree- 
ment within the uncertainties. IMarino et al.1 have dis- 
cussed these abundances at length, so we shall not con- 
si der them further. 

IMarino et al.l (|2011bf) detected an enhancement by 
0.10 dex in the [Ca/Fe] ratio in the r + s group of stars 
in M22, which we recover in our data. Other neighbor- 
ing elements not included in that study also exhibit very 
slight differences in our data. To quantify these differ- 
ences, we apply a line-by-line differential analysis, which 
is largely insensitive to uncertainties in the log(g/) val- 
ues and star-to-star systematic effects in the abundance 



Empirical s-process Abundance Distribution 
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TABLE 5 

Mean Abundances in the Three r + s Stars 



Species 


Z 




1-27 










1-53 








1-80 








(log e) 


{[X/bc\) 


YV 


cr 


o> 


<loge> 


([X/Fej 


) jv 


cr 


°> 


<l°ge} 


{[X/bc\) 


_/V 




cr 


cr M 


Na I 


11 


4.75 


0.39 


3 


0.16 


0.092 


4.70 


0.27 


■1 


0.08 


0.039 


5.11 


0.62 


4 





.06 


0.029 


Mg I 


12 


6.13 


0.41 


3 


0.26 


0.150 


6.28 


0.49 


3 


0.15 


0.085 


6.25 


0.40 


3 





.18 


0.105 


Al I 


13 


5.02 


0.45 


1 


0.23 


0.115 


5.09 


0.45 


■1 


0.13 


0.063 


5.57 


0.87 


1 





.10 


0.048 


Si I 


1 1 


6.02 


0.39 


3 


0.09 


0.053 


5.88 


0.18 


3 


0.21 


0.123 


5.94 


0.18 


3 





.19 


0.109 


K I 


19 


4.14 


0.99 


1 


0.11 


0.110 


4.09 


0.87 


1 


0.11 


0.110 


4.12 


0.84 


1 





.11 


0.110 


Ca I 


20 


5.00 


0.54 


8 


0.10 


0.035 


4.99 


0.46 


8 


0.08 


0.028 


5.04 


0.45 


8 





.11 


0.038 


Sc II 


21 


1.73 


0.38 


5 


0.29 


0.132 


1.60 


0.24 


5 


0.28 


0.124 


1.66 


0.26 


5 





.33 


0.147 


Ti I 


22 


3.05 


-0.02 


9 


0.20 


0.066 


3.26 


0.12 


9 


0.14 


0.048 


3.26 


0.06 


9 


0. 


.12 


0.040 


Ti II 


22 


3.53 


0.39 


9 


0.14 


0.048 


3.65 


0.49 


9 


0.14 


0.046 


3.66 


0.46 


9 





.13 


0.042 


V I 


23 


2.06 


0.02 


5 


0.15 


0.066 


1.98 


—0.14 


5 


0.11 


0.051 


2.02 


—0.16 


5 


0. 


.16 


0.072 


Cr I 


24 


3.54 


—0.22 


6 


0.17 


0.067 


3.67 


—0.16 


6 


0.06 


0.026 


3.71 


—0.18 


6 





.07 


0.030 


Cr II 


24 


3.99 


0.15 


1 


0.11 


0.110 


4.05 


0.20 


1 


0.11 


0.110 


4.26 


0.37 


1 





.11 


0.110 


Mn I 


25 


3.10 


-0.45 


1 


0.16 


0.080 


3.08 


-0.54 


1 


0.10 


0.048 


3.12 


-0.56 


1 


0. 


.13 


0.063 


Fe I 


26 


5.62 


— 1.88 


69 


0.17 


0.020 


5.68 


—1.82 


69 


0.12 


0.015 


5.75 


— 1.75 


69 





.13 


0.016 


Fe II 


26 


5.70 


-1.80 


7 


0.12 


0.046 


5.71 


—1.79 


7 


0.13 


0.048 


5.75 


-1.75 


7 





.14 


0.054 


Co I 


27 


3.17 


0.06 


1 


0.08 


0.039 


3.06 


—0.12 


1 


0.13 


0.064 


3.15 


-0.09 


1 





.09 


0.045 


Ni I 


28 


4.24 


—0.10 


10 


0.10 


0.032 


4.23 


—0.17 


10 


0.08 


0.027 


4.31 


—0.16 


10 





.08 


0.025 


Cu I 


29 


1.95 


—0.36 


2 


0.08 


0.057 


1.87 


—0.51 


2 


0.08 


0.057 


1.84 


—0.60 


2 





.08 


0.057 


Zn I 


30 


2.89 


0.21 


2 


0.08 


0.057 


3.03 


0.28 


2 


0.11 


0.080 


3.08 


0.27 


2 





.10 


0.070 


Rb I 


37 


<1.30 


<0.66 


1 






<1.20 


<0.49 


1 






<1.35 


<0.58 


1 








Sr I 


38 


1.08 


0.09 


1 


0.11 


0.110 


1.08 


0.02 


1 


0.11 


0.110 


0.97 


-0.15 


1 





.11 


0.110 


Y II 


39 


0.77 


0.36 


7 


0.08 


0.031 


0.69 


0.27 


7 


0.05 


0.017 


0.72 


0.26 


7 





.09 


0.035 


Zr I 


40 


1.07 


0.37 


2 


0.08 


0.057 


1.18 


0.40 


2 


0.11 


0.075 


1.10 


0.27 


2 





.08 


0.057 


Zr II 


10 


1.47 


0.69 


3 


0.28 


0.160 


1.31 


0.52 


3 


0.07 


0.041 


1.43 


0.60 


3 





.12 


0.068 


Mo I 


12 


0.32 


0.32 


3 


0.45 


0.259 


0.37 


0.30 


3 


0.55 


0.388 


0.18 


0.05 


2 





.23 


0.130 


Ru I 


4-1 


0.15 


0.28 


2 


0.08 


0.057 


0.00 


0.06 


1 


0.14 


0.100 


0.50 


0.50 


1 





.14 


0.100 


Rh I 


15 


—0.45 


0.37 


1 


0.11 


0.110 


—0.55 


0.20 


1 


0.20 


0.200 


—0.80 


—0.11 


1 





.11 


0.110 


Ba II 


56 


1.20 


0.82 


2 


0.08 


0.057 


1.13 


0.74 


2 


0.08 


0.057 


1.09 


0.66 


2 





.08 


0.057 


La II 


57 


—0.11 


0.59 


14 


0.14 


0.038 


—0.13 


0.56 


1 1 


0.07 


0.020 


—0.21 


0.44 


14 





.08 


0.021 


Ce II 


58 


0.25 


0.48 


16 


0.13 


0.032 


0.43 


0.64 


16 


0.12 


0.029 


0.21 


0.38 


16 





.11 


0.028 


Pr II 


59 


—0.60 


0.49 


1 


0.06 


0.029 


—0.61 


0.47 


4 


0.06 


0.030 


—0.72 


0.31 


1 





.06 


0.030 


Nd II 


60 


0.22 


0.61 


24 


0.10 


0.020 


0.19 


0.56 


21 


0.09 


0.019 


0.06 


0.39 


24 





.09 


0.019 


Sm II 


62 


—0.38 


0.46 


9 


0.12 


0.039 


—0.33 


0.50 


9 


0.09 


0.029 


—0.38 


0.41 


9 





.09 


0.029 


Eu II 


63 


-1.05 


0.23 


3 


0.26 


0.149 


-0.91 


0.36 


3 


0.13 


0.073 


-0.91 


0.32 


3 





.11 


0.065 


Gd II 


64 


—0.25 


0.48 


3 


0.11 


0.065 


—0.15 


0.57 


3 


0.12 


0.070 


—0.32 


0.37 


2 





.22 


0.155 


Tb II 


65 













-1.28 


0.21 


1 


0.11 


0.110 


-1.30 


0.15 


1 





.11 


0.110 


Dy II 


66 


-0.25 


0.46 


3 


0.14 


0.079 


-0.08 


0.61 


4 


0.13 


0.067 


-0.16 


0.49 


3 





.06 


0.035 


Ho II 


67 


-1.15 


0.17 


1 


0.20 


0.200 


-1.18 


0.13 


1 


0.20 


0.200 


-1.20 


0.07 


1 





.20 


0.200 


Er II 


68 


-0.41 


0.47 


2 


0.16 


0.110 


-0.43 


0.45 


2 


0.19 


0.135 


-0.25 


0.59 


2 





.13 


0.095 


Tm II 


69 













-1.50 


0.19 


1 


0.20 


0.200 















Yb II 


70 


-0.76 


0.12 


1 


0.11 


0.110 


-0.56 


0.31 


1 


0.11 


0.110 


-0.52 


0.31 


1 





.11 


0.110 


Hf II 


72 


-0.56 


0.39 


1 


0.11 


0.110 


-0.55 


0.39 


1 


0.11 


0.110 


-0.60 


0.30 


1 





.11 


0.110 


Ir I 


77 





































Pb I 


82 


0.77 


0.61 


1 


0.11 


0.110 


0.97 


0.74 


1 


0.11 


0.110 


0.63 


0.34 


1 





.11 


0.110 


Th II 


90 





































Note. 


— Quoted uncertainties represent internal uncertainties only. {[Fe/H]) 


is listed 


in the ([X/Fe]} column for Fe I and Fe II. 





emcnts heavier than Fe may not be surprising, since 
s-process nucleosynthesis produces heavy nuclei from 
successive neutron capture on Fe-group seeds. Variations 
in the lighter Fe-group elements are more surprising. We 
return to this issue in Section [6] 

4.2. The Neutron- Capture Abundance Patterns in M22 

Figure [3] illustrates the abundance patterns for the 
Z > 38 elements in each of the 6 stars observed 
in M22. The abundance pattern of the metal- 
poor ([Fe/H] = —2.1) r-process-rich standard star 
BD+17 3248 ([Eu/Fc] = +0.9) is shown for compari- 
son. The three stars selected from the s-poor group of 
iMarino et" al. (2011b) — our "r-only" group — share a sim- 
ilar abundance pattern with each other and BD+17 3248. 
The three stars from the IMarino et al.l s-rich group — our 
"r + s group" — share a similar abundance pattern with 
each other that clearly differs from BD+17 3248 for the 
lighter n-capture elements and Pb. Figure [3] demon- 



strates that it is appropriate to average together the 
abundances of the 3 stars in each of these two groups 
to reduce random uncertainties in the abundances, par- 
ticularly for the abundances derived from small numbers 
of lines. The average abundance patterns for the r-only 
and r + s groups are shown in Figure 0J The derived 
mean [Fe/H] for the 3 stars in each group is the same, 
so the relative vertical scaling of the abundances in Fig- 
ure [4] is not affected by the bulk metal content of these 
two groups. 

In the 7'-only group, the abundance pattern for Ba and 
the heavier elements (Z > 56) generally conforms to that 
of BD+17 3248. When normalized to Eu (Z = 63), the 
Ba, Ce, and Nd [Z = 56, 58, and 60, respectively) abun- 
dances in the r-only group appear slightly enhanced rel- 
ative to BD+17 3248. Furthermore, in the r-only group 
several of the odd-Z elements in the rare earth domain 
(Tb, Ho, and Tm — elements 65, 67, and 69, respectively) 
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Fig. 3. — Logarithmic abundances for Z > 38 elements in the 3 r-only stars (blue crosses, left panels) and the 3 r + s stars (red 
circles, right panels ) as a function of atomic number. The gray line illustrates the abundances in the r-process standard star BD+17 3248 
ICowaji et al.lr2002l . [20051; ISneden et a,l.|[2009l : IRoederer et al-lfeoTOcTl . Pb h as not been detected in BD+17 3248, so we instead show the 
predicted Pb/Eu ratio based on the average Pb/Eu observed in Figure 3 of Roederer et al. (2010b). The BD+17 3248 abundance pattern 
has been normalized to the Eu abundance in each star. 



plus the even-Z element Yb (Z = 70) lie 0.2-0.4 dex 
below the BD+17 3248 abundances. This is not sur- 
prising given that the r-process enrichment in M22 is 
less extreme than that seen in BD+17 3248 or other r- 
rich standard stars, and variations in the physical con- 
ditions at the time of the nucle osynthesis may be re- 
sponsible (Ro ederer et al.l [2010b| ). The abundances of 
the lighter elements Sr-Rh (38 < Z < 45) are known 
to vary widely among metal-poor star s that show no 
eviden ce of s-process enrichment (e.g., IRoederer et al.l 
l2010bl and references therein). Based on the empiri- 
cal correlation betwee n [Eu/Yl and [Eu/F e l identified by 
Barkl em et al.1 (120051) .lOtsuki et all (|2006Q . Montes etall 
(12007ft . and lRoederer et al.l (|20ldbf ) . we would expect the 
Sr-Rh elements in the M22 r-only group to be more 
abundant than those in BD+17 3248 when normalized 
to Eu, which is indeed the case. These elements may be 
produced by primary nucleosynthetic mechanisms in ad- 
dition to the r-process (e.g., charged-particle reactions 
in the expanding neutrino winds of core collapse SNe; 
iWooslev k. Hoffman! 11992ft and so could be expected to 
vary. 

In the r + s group, all heavy elements except Mo 
(Z = 42), Ru (Z = 44), Eu, Ho, and Tm are en- 
hanced relative to the r-only group. These differences 
are most pronounced among the lightest n-capture ele- 
ments (Sr, Y, and Zr), the light and heavy ends of the 
rare earth domain (Ba-Nd and Yb-Hf), and Pb. This 
is not surprising, given that a significant fraction of the 
S.S. abundance of each of these elements is attributed to 



the s-process. In contrast, the S.S. abundances of ele- 
ments in the middle of the rare earth domain are mostly 
attributed to the 7'-process. 

Low-metallicity AGBs produce substantial overabun- 
dances of Pb relative to the Fe-group s-process seeds 
and all elements intermediate betw een Fe and Pb (e.g., 
iClavtonl Il98l iGallino et all 11998ft . As the metallicity 
of the s-process environment increases above [Fe/H] ~ 
— 1.0 , the Pb overabundances decrease ()Travaglio et al.l 
120011 ). IRoederer et~aT1 (|2010bD have shown that [Pb/Eu] 
ratios can be an effective diagnostic to identify low metal- 
licity stars that lack detectable contributions from the 
s-process. It is clear from Figures [TJ [3j and [4] that the 
Pb abundance is moderately enhanced in the r + s group 
of stars relative to the r-only group. As shown in Fig- 
ure [La/Eu] and [Pb/Eu] in M5, M 13 M15, M92 
and N GC 6752 (lYong et al ll2fM [2008313; ISobeck et al.l 
l2lmt IRoederer fc Snedenl 1201 lb are the same as that 
for field stars of the same metallicity. These ratios are 
low and suggest no contribution from s-process material. 
The M22 r-only group is normal for other metal-poor 
GCs in this regard. [Pb/Eu] is moderately enhanced in 
the M22 r + s group, and this increase relative to the r- 
only group (a difference of +0.85 dex) is notably higher 
than other [X/Eu] rat ios (< +0.55 dex). This fu rther 
confirms the results of iMarino et al.l (|2009l l2011bft that 
the r + s (or s-rich) group in M22 contains a moderate 
amount of s-process material. 
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Fig. 4. — Top panel: The mean logarithmic abundances for the 3 r-only stars (blue crosses) and the 3 r + s stars (red circles) as a 
function of atomic numb e r. The gray line and small gray squares i llustrate the abundances in the r-process standard star BD+17 3248 
llCowan et al.ll2003 , [20051; ISneden et al.ll2009l ; IRoederer et alj|2uTDcD . Pb h as not been detected in BD+17 3248, so we instead show the 
predicted Pb/Eu ratio based on the average Pb/Eu observed in Figure 3 of Rocdercr et al. (2010b). The BD+17 3248 abundance pattern 
has been normalized to the Eu abundance. Bottom panel: The differences in these mean abundances. The dotted line indicates zero 
difference. 



4.3. The Age of M22 Calculated from Radioactive 
232 Th Decay 

The radioactive isotope 232 Th can only be produced 
in r-process nucleosynthesis. It can be used in con- 
junction with other stable elements produced in the 
same events to yield an age for the r-process material 
in M22. This can be done in a relative sense (e.g., 
comparing the Th/Eu ratio in several GCs) or an ab- 
solute sense if the initial production ratio of Th/Eu 
is known from theory. We u se the production ratio 
predicted by the simulations of iKratz et al.l (|2007t ) and 
the derived loge(Th/Eu) ratio in the 3 r-only stars in 
M22 (-0.60 ± 0.085) to calculate an absolute age of 
12.4 ± 4.0 Gyr. Recall we could not measure Th in 
the r + s group due to blending features. This assumes 
no uncertainty in the initial production ratio, which 
likely translates to an uncertainty of several Gyr (e.g . , 
iFrebel et al.ll2007t IKratz et al.|[2007l: ILudwig et anl2010l ) . 
This age estimate is consistent with the relatively old 
age derived fro m isochrone fitting to the M22 main se- 
quence turnoff (jMarm-Franch et ail I2009H ■ the ages of 
other metal-poor GCs derived from their Th/Eu ratios 
( Sne den et al.l 120001 [Johnson fc Boltd I2001L lYong et al.l 
2008bl lLai et al.H201in 7 and halo field stars of similar low 



metallicity (e.g.. IRoederer et aL1l2009|) . While the useful- 
ness of this measurement is limited by observational un- 
certainties and systematic effects, the general agreement 
is reassuring. 

5. COMPARISON TO OTHER COMPLEX 
METAL-POOR GCS 

As discussed in IMarino et al.l (|2009t 1201 lB . evidence 
for multiple stellar populations in M22 includes the fol- 
lowing: (1) the SGB shows two distinct sequences, (2) 
there is a metallicity offset between the two groups, 
(3) each group independently exhibits the O-Na and 
C-N anticorrelations, and (4) there are clearly distinct 
n-capture abundance patterns in the two groups. It is 
difficult to envision an unambiguous evolutionary pic- 
ture for M22 that accounts for the entire body of ob- 
servations. Here, we illuminate this issue by comparing 
M22 with other GCs that show similar complexity, like 
NGC 1851, and simpler GCs, like M4 and M50 

The GCs M4 and M5 are a frequently studied pair 

9 M22 is among the more massive Milky Way GCs 
(4.0 X 10 5 Mq, assuming M/L v = 2 Mq/Lq), and the present- 
day mass of M22 is also similar to that of M4, M5, and NGC 1851 
(1.2 X 10 5 Mq, 5.4 X 10 5 Mq, and 3.4 X 10^ Mq, respectively). 
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Fig. 5. — [La/Eu] and [Pb/Eu] ratios as a function of [Fe/H]. Only GCs where Pb has been measured have been included. The dotted 
lines indicate the solar ratio. A typical uncer tainty is shown. The blue long-dashed lines indicate the approximate yields of a 5 Mq AGB 
star at [Fe/H] = —2.3 (Roederer ct al. 201013) ■ The red short-dashed lines indicate the approximate mean s of metal-poor field stars whose 
[Pb/Eu] ratios are cons istent with having been enriched by r-proce ss material o nly (Roederer ct al. 2010b). GC abund ances are referenced 
as foll ows: M4 and M5. IYone eFaLI |2008al|B); M13 a nd NG C 6752, lYong et""aLl H20091 ; M15. ISobeck et"aTH20Tiri ; M92. IRoederer fc Snedenl 
(2011); M22, this study; field stars. IRoederer eF al. (2010b). All abundances have been normalized to the scale used in the present study. 



of clusters that are not physically related to one an- 
other. Both are more metal-rich than M22 ([Fe/H] = 
— 1.2 and —1.3), and previous work has revealed that 
M4 con tains moderate s - proces s enrichment relati ve 
to M5 divans et all H99l I200H: lYong et all [2008allbl : 
iMarino et al.ll2008f) . The heavy element abundances in 
M5 are similar to t he scaled S.S. r -process residuals 
(|Yong et all l2008alTbl: lLai et aHl2lJTTI ). and the low Pb 
abundance (jYong et al.ll2008ai r suggests that these ele- 
ments were produced by r-process nucleosynthesis with- 
out the need to invoke contributions f rom the s-process 
(jRoederer et al.ll2010rjrRoedererll2011l) . We subtract the 
heavy element abundances in M5 from those in M4 (cf. 



lYong et all l2008bl ) to estimate the s-process contribu- 
tion to M4. As shown in Figure [6j these differences 
are remarkably similar to the differences observed be- 
tween the r + s and r-only groups in M22. There is a 
gradual increase in the s-process content of Co through 
Zn (27 < Z < 30), a moderate s-process contribu- 
tion with some element-to-element scatter for Rb-Rh 
(37 < Z < 45), a gradual decrease from Ba to Gd 
(56 < Z < 64), and a gradual increase from Yb to Pb 
(70 < Z < 82) 03 There is no a priori reason to expect 

10 Neither llvans et al.l 120011) nor lYong etaT] (I2008bl) found dif- 
ferences in [Ca/Fe] between M4 and M5. Those studies did not 
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TABLE 6 

Mean Abundances in the r and r + s Groups 



r-only 



r + ; 



Species 


z 




°> 




°> 


Na I 


11 


0.08 


0.018 


0.49 


0.023 


Mg I 


12 


0.39 


0.064 


0.45 


0.060 


Al i 


13 


0.34 


0.034 


0.69 


0.036 


Si I 


14 


0.14 


0.064 


0.33 


0.044 


K I 


19 


0.80 


0.064 


0.90 


0.064 


Ca I 


20 


0.40 


0.023 


0.48 


0.019 


Sc II 


21 


0.21 


0.073 


0.30 


0.077 


Ti I 


22 


0.06 


0.015 


0.07 


0.028 


Ti II 


22 


0.47 


0.020 


0.45 


0.026 


V I 


23 


-0.14 


0.026 


-0.10 


0.035 


Cr I 


24 


-0.14 


0.020 


-0.17 


0.019 


Cr II 


24 


0.14 


0.064 


0.24 


0.064 


Mn I 


25 


-0.50 


0.034 


-0.53 


0.034 


Fe I 


26 


-1.81 


0.008 


-1.81 


0.009 


Fe II 


26 


-1.78 


0.022 


-1.78 


0.028 


Co I 


27 


-0.15 


0.033 


-0.03 


0.027 


Ni I 


28 


-0.20 


0.020 


-0.15 


0.016 


Cu I 


29 


-0.66 


0.033 


-0.49 


0.033 


Zn I 


30 


0.04 


0.032 


0.24 


0.039 


Rb I 


37 


<0.20 




<0.49 




Sr I 


38 


-0.53 


0.064 


-0.01 


0.064 


Y II 


39 


-0.24 


0.013 


0.29 


0.014 


Zr I 


40 


-0.08 


0.036 


0.34 


0.035 


Zr II 


40 


0.22 


0.030 


0.55 


0.034 


Mo I 


42 


0.11 


0.078 


0.12 


0.111 


Ru I 


44 


0.13 


0.053 


0.28 


0.044 


Rh I 


15 


-0.30 


0.078 


0.14 


0.072 


Ba II 


56 


0.18 


0.047 


0.74 


0.033 


La II 


57 


0.08 


0.016 


0.51 


0.013 


Ce II 


58 


-0.07 


0.009 


0.49 


0.017 


Pr II 


59 


0.12 


0.016 


0.42 


0.017 


Nd II 


60 


0.15 


0.009 


0.52 


0.011 


Sm II 


62 


0.27 


0.013 


0.46 


0.018 


Eu II 


63 


0.35 


0.056 


0.32 


0.046 


Gd II 


64 


0.33 


n 04*3 
u.u^o 


0.51 


0.046 


Tb II 


65 


-0.05 


0.064 


0.18 


0.078 


Dy II 


66 


0.42 


0.036 


0.51 


0.029 


Ho II 


67 


0.12 


0.115 


0.12 


0.115 


Er II 


68 


0.38 


0.033 


0.52 


0.063 


Tm II 


69 


0.09 


0.045 


0.19 


0.200 


Yb II 


70 


-0.11 


0.115 


0.25 


0.064 


Hf II 


72 


-0.11 


0.064 


0.36 


0.064 


Ir I 


77 


0.45 


0.064 






Pb I 


82 


-0.26 


0.064 


0.56 


0.064 


Th II 


90 


0.21 


0.064 






Note. 


— The means 


represent 


weighted means from 



the 3 stars in each group, and the stated uncertainties 
represent internal uncertainties only. ([Fe/H]} is listed 
in the ([X/Fe]} column for Fe I and Fe II. 



such similarity. Figure [5] implies that the heavy elements 
in M5 and the M22 r-only group were produced by simi- 
lar nucleosynthesis mechanisms, and the heavy elements 
in M4 and the M22 r + s group were produced by anothcr 
similar set of nucleosynthesis mechanisms. 

The heavy eleme nts in NGC 1851 resemble the pattern 
obser ved in M22 (lYong fe Grundahl 120081: lYong et al ' 
2009; iCarretta et al.l 120101 1201 ll) . and iCarretta et~a! 



raised the possibility that NGC 1851 may have formed 
through the merger of two proto-clusters in a now dis- 
solved dwarf galaxy. The Eu abundance within each of 
M22 and NGC 1851 is constant, but moderate enhance- 
ments are o bserved in Zr, B a, La, and Ce in some stars of 
both GCs. ICarretta et al.l report a small but detectable 

examine K, and the rest of the abundance ratios from Ca-Mn in 
M4 and M5 were found to be identical. 
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Fig. 6. — Top panel: Differences between the mean abundances 
in GCs M4 a nd M5. The abundances are taken from lYong et al.l 
H2008allbl) and llvans eTaTI SiWH 1200111 . Bottom panel: Differences 
between the mean abundances in the r + s and r groups in M22. 
In both panels, dotted lines indicate zero difference. 



TABLE 7 

Mean Line-by-line Differentials for 
K-Zn in the r and r + s Groups 



Species 


JV 


<A[X/Fe]) a 


(7 


°> 


K I 


1 


+0.100 


0.090 


0.090 


Ca I 


8 


+0.066 


0.063 


0.022 


Sc II 


5 


+0.087 


0.048 


0.022 


Ti I 


9 


-0.007 


0.087 


0.029 


Ti II 


9 


-0.003 


0.062 


0.021 


Ti I+II 


18 


-0.005 


0.073 


0.017 


V I 


5 


+0.057 


0.060 


0.027 


Cr I 


6 


-0.016 


0.073 


0.030 


Cr II 


1 


+0.105 


0.110 


0.110 


Cr I+II 


7 


-0.008 


0.085 


0.032 


Mn I 


4 


-0.064 


0.078 


0.039 


Co I 


4 


+0.094 


0.065 


0.033 


Ni I 


10 


+0.029 


0.073 


0.023 


Cu I 


2 


+0.169 


0.088 


0.062 


Zn I 


2 


+0.189 


0.052 


0.036 



In the sense of ([X/Fe] r+S ) - <[X/Fe] r ) 



spread in Fe and Ca in NGC 1851. Unlike M22, these 
two elements are strongly correlated, which implies the 
more metal- rich group is not enhanced in [Ca/ Fe] relative 
to the metal-poor group. Note, however, that iLee et all 
(2009) suggest much larger [Ca/H] variations are present 
in NGC 1851 Lik e M22, NGC 1851 has a split SGB 
(jMilone et al.|[200l ). which may be explained by either 
an age difference of ~ 1 Gyr or a differenc e in the over- 
all CNO with a negligible age difference (jCassisi et al.l 
120081 I Ventura et aITl2009( ). Examination of the radial 
distributions of different SGB populations gives conflict- 
ing results for NGC 1851, and radial distributions for the 
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two groups of stars in M22 have not been investigated. 

The M22 chemistry does not exclude the possibil- 
ity that it formed through a merger of two separate 
groups similar to M4 and M5 (at lower metallicity). 
The s-process abundances in the two M22 groups are 
sufficiently distinct that these two groups would be re- 
garded as completely separate populations if not ob- 
served together in the same GC. Similar metallicities and 
r-process abundances might be expected if the groups 
formed in close proximity in a now-dissolved dwarf 
galaxy. 

On the other hand, M22 shares several characteris- 
tics with the metal-poor populations in uj Cen, which 
is more difficult to interpret as having formed via merg- 
ing of several clusters. Based on current self-enrichment 
models, a possible way to account for the M22 chemistry 
is through fine-tuning of the times of accumulation of 
the material from which successive generations form. In 
this scenario M22 does not evolve as an isolated system, 
and extern al gas flows can cont ribute to the enrichment 
processes dMarino et al.l l2011b| ). A similar mec h anism 
has been recently suggested by iD'Antona et al.l (|201lD 
to explain the O-Na antic orrelation pattern in the more 
complex case of u Cen ([Johnson fc Pilachowskil 120101 : 
iMarino et al.ll20~lla|) . Further exploration is beyond the 
scope of the present work. 

6. THE SOURCE OF THE S-PROCESS MATERIAL 

To summarize the results of the previous sections, the 
heavy elements in the M22 r-only group can be explained 
by nucleosynthesis mechanisms associated with core col- 
lapse SNe. The r + s group contains a moderate amount 
of material produced by s-process nucleosynthesis. Pre- 
vious studies have shown that the r+s group has a higher 
mean metallicity than the r-only group, but the ratio of 
r-process material to Fe-group material is roughly equal 
in the two groups. In this section we investigate possible 
nucleosynthetic sources for the s-process material in the 
r + s group. 

We subtract the r-process contribution (i.e., the abun- 
dance in the r-only group) to each of La and Pb in 
the r + s group to derive the intrinsic [Pb/La] s ratio. 
We perform a similar calculation to estimate the intrin- 
sic s-process rati os in M4 by subtra cting the M5 abun- 
dances using the lYong et all l|2008af ) abundances. This 
yields [Pb/La] s = +0.18 ± 0.09 in M22 and [Pb/La] s = 
-0.01 ± 0.08 in M4. Similarly, we derive the indice£3 
[hs/ls} s = -0.01 and -0.50 and [Pb/hs] = +0.29 and 
+0.28 for M22 and M4, respectively. (Uncertainties on 
each of these quantities are likely 0.10-0.15 dex.) These 
ratios and indices are useful since they are insensitive to 
the dredge-up efficiency or the dilution of AGB products 
in the stars currently observed. We infer that the AGBs 
providing the s-process enrichment in M4 and the r + s 
group in M22 were similar but not identical. 

Models of s-process nucleosynthesis indicate that Pb is 
a sensitiv e probe of the stellar mass, metallicity, and neu- 
tron flux. iGoriely fc Mowlavil (|2000l ) present yields for a 
model representative of 1.5 < M < 3.0 M , [Fe/H] = 

11 As defined by, e.g., IBisterzo etTal] 112010 ), the 
ratios of light (Is) and heavy (/is) s-process yields 
are [Is/Fe] = i([Y/Fe] + [Zr/Fe]) and [hs/Fe] = 
|([La/Fe] + [Nd/Fc] + [Sm/Fe]). Also, [hs/ls] = [hs/Fe] - [Is/Fe]. 



— 1.25 AG B stars, and [Pb / La] c an be estimated from 
Figure 3 of IGoriely fc Siessl (120011) for th eir 3 zero- 
metallicity AGB model. ICristallo et al.l (|2009l ) present 
yields for 2 M@ AGB mod els at [Fe/H] = -1.2 and 
—2.2. IBisterzo et al.l (|201C0 present a set of yields for 
several masses (M = 1.3, 1.4, 1.5, and 2.0 M ), metal- 
licities (-3.6 < [Fe/H] < -1.0), and 13 C pocket ef- 
ficiencies. [Pb/La] predictions can also be calculated 
for limited combinations of masses, metallicities, and 
13 C pocket efficiencies fro m the AGB yields presented 
in iRoederer et alj ()2010bD . These predictions rely on 
similar atomic data, stellar models, assumptions about 
branching points, etc., and so are not entirely indepen- 
dent. 

When compared with these yields, the M4 and M22 
s-process heavy element ratios and indices point to a 
common theme: low mass AGB stars (M < 3 M ) can- 
not reproduce the observed values unless the standard 
13 C pocket efficiency is reduced by factors of 30-150. Pb 
is enhanced in both M4 and the r + s group in M22 rel- 
ative to the lighter n-capture elements and Fe, but it is 
not nearly as enhanced as observed in metal-poor stars 
extrinsically enriched in s-process elements by an AGB 
binary companion. AGBs with M ~ 4.5-6.0 M (those 
which may not form a 13 C pocket and hence will not acti- 
vate the 13 C(a, n) 16 Q neutron source) can produce lower 
[Pb/La] ratios (jRoederer et al.ll2010bl) . For comparison, 
predictions for the 5 M AGB mo dels at [Fe/H] = -2.3 
are shown in Figured] Figures in IBisterzo et al.l (|2010f ) 
present the [hs/ls] and [Pb//is] indices for a limited num- 
ber of 3 and 5 M AGB models. Their predictions for 
the appropriate (low) 13 C pocket efficiency in a 5 M 
AGB are a near-perfect match to the s-process ratios 
in each of M4 and M22 at their respective metallicities. 
This result is encouraging. 

The 22 Ne neutron source, which activates at higher 
temperatures than the 13 C neutron source, does not 
play a dominant role in AGB stars with M < 3- 
4 M Q . In AGB stars with M = 5-8 M , the tem- 
perature at the base of the thermal pulse is higher, 
and the 22 Ne(a,n ) 25 Mg reaction can occur there (e.g., 
iBusso et al.ll200ll ). In principle this could also account 
for the s-process neutron captures that produ c e smal l 
amo unts of Co-Zn , as ob served; see lYong et all (|2008bn 
and iKarakas et all (|2009f ) for further discussion. 

Models of the weak component of the s-process 
have traditionally been set in ~ 25 M stars (e.g., 
iRaiteri et al.lll993T ) that activate the 22 Ne neutron source 
during core He-burning and shell C-burning stages, since 
models of less massive stars suggest that subsequent 
burning stages will destroy any s-process material cre- 
ated. Models that include rotationally-induced mixing 
can increase the neutron flux by mixing 14 N (which is 
converted to 22 Ne) into the relevant regions, possibly pro- 
ducing nuclei as heavy as 208 Pb (jPignatari et al.l l2008). 
Yet the enhanced s-process abundances observed in the 
M22 r + s group cannot be due to the operation of the 
weak s-process in massive stars. There is no reason to ex- 
pect that the the SNe that enriched the metal-rich r + s 
group in M22 host the weak s-process and those that 
enriched the metal-poor r-only group did not. 

The minority neutron-rich Mg isotopes 25 Mg 
and 26 Mg may be produced (among other proton- 
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TABLE 8 

T- AND S-PROCESS PERCENTAGES IN THE M22 r + S GROUP 
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0.92 


30.2 


69.8 


Y 


39 


0.151 


0.0447 
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0.72 


0.93 
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62.0 
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1.35 


1.02 


1.08 


46.8 


53.2 


Mo 


12 


0.0479 


0.0468 


0.00109 


0.22 


0.21 


-1.42 


97.7 


2.3 


Ru 


11 


0.0513 


0.0363 


0.0150 


0.25 


0.10 


-0.28 


70.8 


29.2 


Rh 


15 


0.00759 


0.00275 


0.00483 


-0.58 


-1.02 


-0.78 


36.3 


63.7 


Ba 


56 


0.398 


0.110 


0.288 


1.14 


0.58 


1.00 


27.5 


72.5 


La 


57 


0.0195 


0.00724 


0.0122 


-0.17 


-0.60 


-0.37 


37.2 


62.8 


Ce 


58 
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0.0155 


0.0408 


0.29 


-0.27 


0.15 


27.5 


72.5 


Pr 


59 


0.00661 


0.00331 


0.00330 


-0.64 


-0.94 


-0.94 


50.1 


49.9 


Nd 


60 


0.0417 


0.0178 


0.0239 


0.16 


-0.21 


-0.08 


42.7 


57.3 


Sm 


62 


0.0126 


0.00813 


0.00446 


-0.36 


-0.55 


-0.81 


64.6 


35.4 


Eu 


63 


0.00331 


0.00355 


-0.000237 c 


-0.94 


-0.91 




100. 
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Gd 


64 
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Dy 


66 


0.0195 


0.0158 


0.00365 


-0.17 


-0.26 


-0.90 


81.3 


18.7 


Ho 


67 


0.00191 


0.00191 


0.00 


-1.18 


-1.18 




100. 


0.0 


Er 


68 


0.0132 


0.00955 


0.00363 


-0.34 


-0.48 


-0.90 


72.4 


27.6 


Tm 


69 


0.000933 


0.000741 


0.000192 


-1.49 


-1.59 


-2.18 


79.4 


20.6 


Yb 


70 


0.00708 


0.00309 


0.00399 


-0.61 


-0.97 


-0.86 


43.7 


56.3 


Hf 


72 


0.00776 


0.00263 


0.00513 


-0.57 


-1.04 


-0.75 


33.9 


66.1 


Ir 


77 




0.0324 






0.05 








Pb 


82 


0.191 


0.0288 


0.162 


0.82 


0.00 


0.75 


15.1 


84.9 


Th 


90 




0.000891 
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+ 11.2 
-9.8 
+7.7 
-7.1 
+11.1 
-9.5 
+ 11.2 
-9.5 
+ 10.4 
-7.7 
+ 1.9 
-1.8 
+6.8 
-5.7 
+7.4 
-6.4 
+53.3 
-2.3 
+ 17.7 
-14.2 
+ 15.0 
-10.6 
+5.6 
-4.6 
+2.6 
-2.4 
+ 1.7 
-1.6 
+4.0 
-3.7 
+2.0 
-1.9 
+4.8 
-4.4 
+28.4 
-0.0 
+ 15.0 
-12.2 
+22.8 
-16.4 
+ 13.1 
-11.3 
+69.8 
-0.0 
+ 17.9 
-14.4 
+60.2 
-20.6 
+22.3 
-14.7 
+ 11.6 



+5.2 
-3.9 



a The r component implicitly includes contributions from all other processes that may have enriched the 
stars in M22 prior to the epoch of s-process enrichment, e.g., charged-particle reactions, etc. 
b loge = logA r + 1.54 

c Indicates mild destruction of Eu by the s-process (not statistically significant) 



and a-capture channels) by the reaction sequence 
22 Ne(a,n) 25 Mg(n,7) 26 Mg, which acts as both a neu- 
tron source and poison. Preliminary measurements 
of ( 25 Mg+ 26 Mg)/ 24 Mg in M4 and M5 indicate that 
the Mg isotopes have sim ilar proportions in the two 
clusters (jYong et al.1 1200 8b). While preliminary, these 
measurements hint that the source affecting the Mg 
isotopic ratios has acted similarly in M4 and M5. Since 
moderate quantities of s-process material are observed 
in M4 and M22 but not M5, it seems unlikely that the 
source of the s-process material modifies the Mg isotopic 
ratios substantially. Unfortunately we cannot assess 
the Mg isotopic ratios from our M22 data, but new 
measurements of these ratios in all three GCs would be 
of great interest. 

At low metallicity 22 Ne also serves as a primary seed 
nucleus from which a chain of n-capture reactions can 
generate a s mall leakage across the Fe-group isotopes 
(jBusso et al]l27MlGalhno et al.ll2006h . We propose that 
the observed variations in the Fe-group ratios and per- 
haps even the overall metallicity (Fe) increase in the r + s 
group could be due to this phenomenon. The s-process 
path passes through stable or long-lived nuclei of K, Ca, 



Ti, V, and Cr, including several nuclei ( 39 K, 42 Ca, 43 Ca, 
44 Ca, 50 Ti, 51 V, 52 Cr) with closed nuclear shells. The 
only stable isotopes of Sc and Mn, 45 Sc and 55 Mn, do 
not have closed nuclear shells, so it is perhaps surprising 
that Sc shows an enhancement while Mn shows a defi- 
ciency in the r + s group. The fact that we observe no 
change in Ti or Cr could be related to the initially larger 
abundances of these even- if elements relative to a small 
s-process contribution. Ca, which could also be expected 
to follow this pattern, may be enhanced because there are 
three Ca isotopes on the s-process path with closed pro- 
ton shells. Obviously detailed calculations are needed to 
test these proposals for the Fe-group variations between 
the two groups in M22. 

If the s-process material in M4 and M22 is produced 
by neutrons from the 22 Ne source, this implies an origin 
different from that of the s-process material in w Cen. 
Smit h~et al.l ([2000) found that the n-capture elements 
in lu Cen are best fit by low-mass (1.5-3.0 M Q ) AGB 
stars where the 13 C neutron source is active. The ob- 
served [Rb/Zr] ratios in u> Cen, which are quite sensi- 
tive to the neutron density and hence the neutron source 
because of s-process branching at 85 Kr, are best fit by 
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Fig. 7. — The fractional component of the r + s group of stars 
originating in the s-process for elements heavier than Fe. 
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low mass (M < 3 M R ) AGB mod els. The [Rb/Zr] 

ratios in M4 derived by lYong et al.l (|2008af) are higher 
than those in uj Cen, and our [Rb/Zr] ratio in the M22 
r + s group is not lower than that in M4, supporting 
our assertion. (Recall that we could only derive up- 
per limits on the R b abundance in M22.) Furthermore, 
ICunha et al.1 {2002) found no evolution in the [Cu/Fe] 
ratio over —2.0 < [Fe/H] < —0.8 in uj Cen, indicat- 
ing that there were no contributions to Cu from AGB 
stars that could produce Cu with neutr o ns fro m the 
22 Ne(a,n) 25 Mg reaction. iD'Antona et al.l (|2011[ ) point 
out that the timescales for establishing the light element 
variations and the s-process enrichcmcnt in uj Cen are 
discrepant, and this issue is not yet resolved. 

These constraints raise an obvious question: if the 
s-process material in M4 and M22 is produced in more 
massive AGB stars, then why is s-process material not 
detected in every clu ster where the light element vari- 
ations are observed? iMarino et al.1 (|2008h showed that 
there might be a weak correlation between [ Ba/Fel 
and [A l/Fe] in M4, a point also investigated by ISmithl 
(2008). Those data also suggest weak correlations be- 
tween [Ba/Fe] and each of [Na/Fe] and [Si/Fe]. For 
the majority of GCs, however, such correlations are no t 
found (e.g., lArmoskv etafl Il994t iD'Orazi et aD I2010T ). 
This supports our conclusion, drawn from the [Pb/Eu] 
ratios, that the heavy elements in most metal-poor GCs 
are produced by r-process nucleosynthesis. Perhaps in 
GCs like M4, the r + s group in M22, or the metal-rich 
group of NGC 1851, material from slightly lower AGB 
masses was allowed to enrich the GC ISM before the 
clusters formed. This might suggest that these particular 
clusters were more massive initially or originated in dwarf 
galaxies whose potentials could more easily retain ejecta 
and sustain extended periods of star formation. This sce- 
nario is appealing because several of the metal-poor clus- 
ters exhibiting s-process enrichment (M22, NGC 1851, 
uj Cen) exhibit at least minimal spreads in Fe and (in 
the case of M22 and NGC 1851) could have been formed 
through mergers. 

In summary, the observed s-process abundance pat- 
terns are not well-fit by low-metallicity models of AGB 



Fig. 8. — Differences in the mean abundances between the r- 
only group and the r + s group as a function of the s-fraction of 
each element in the S.S. The top panel shows the s-fraction as 
calculated from the average yields of the 1. 5 and 3.0 Mq stellar 
models at [Fc/H] = -0.3 of Arlandini et al.l |[T999|), including the 
co ntribution of low - metal licity AGB stars to the S.S. Pb as derived 
by ITrav aglio et al. (200j). The botto m panel shows the s-fr action 
as calculated by the classical method (Simmerer et al. 2004). The 
dotted lines indicate zero difference. 

stars with M < 3 Mq. Higher mass AGBs that acti- 
vate the 22 Ne neutron source may provide a better fit. 
Both stellar groups in M22 exhibit the Na-0 anticorre- 
lation, but the observed lack of a correlation between 
s-process enrichment and Na within the r + s group 
is difficult to understand if these elements are all pro- 
duced by AGB stars of higher masses. We encourage 
more detailed exploration of the possible association be- 
tween the s-process products in these metal-poor GCs 
with intermediate-mass AGB stars. 

7. AN EMPIRICAL S-PROCESS ABUNDANCE 
DISTRIBUTION 

In this section we compare the nature of low-metallicity 
s-process enrichment with the s-process abundance pat- 
tern observed in the S.S. The M22 s-process "residual" 
is derived by subtracting the abundances in the r-only 
group from the abundances in the r + s group. This 
method assumes that the r-process material in both 
groups is identical, as indicated by observations. 

Table |8] lists the abundances and r- and s-process frac- 
tions for the heavy elements in M22. Figure [7] illustrates 
the fraction of each of these elements that originates in 
the s-process in the M22 r + s group. Elements on the 
s-process path with closed neutron shells (Sr-Zr, Ba- 
Nd, and Pb), plus a few others (Rh, Yb, Hf), owe more 
than 50% of their abundance in the r + s group to the 
s-process. Elements in the middle of the rare earth do- 
main (Sm-Tm) and just beyond the first s-process peak 
(Mo, Ru) are still mostly made of r-process material, 
with s-process fractions less than 40% or so. More than 
80% of the Pb in the r + s group originated in the 
s-process, the most of any element studied. Several el- 
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ements, including Mo, Eu, Ho, and Tm, are consistent 
with a pure r-process origin (i.e., show no enhancement 
in the r + s group) within the uncertainties. In analogy 
with S.S. r-residuals derived via the classical approach, 
elements with small s-process fractions have the largest 
s-process fraction uncertainties, and elements with large 
s-process fractions have the smallest uncertainties. 

To compare the s-process fractions in M22 with the 
s-process fractions derived from the stellar model and 
classical approach, Figure |8] displays the elemental abun- 
dance differences between the two M22 groups as a 
function of the s-process fraction in the S.S. Only ele- 
ments produced predominantly by the main and strong 
s-process components in the S.S. are shown (i.e., Z > 56), 
since the yields of these elements should be less sensitive 
to the source of the neutron flux. There is a remark- 
ably clear correlation, which changes littl e when different 
stellar model s-process fr actions (e.g., lArlandini et al.1 
Il999t iBisterzo et al.l 120 10l) or classical method s-process 
fractions (e.g.. iBurris et"aLll200Ct ISimmerer et~al]|2004D 
are used. Note that the s-process fraction of Pb shown 
in Figure [5] accounts for the low-metallicity AGB com- 
ponent a ccording to the Gala ctic chemical evolution 
model of iTravaglio et al.l (1200 If). The ste l lar m odel of 
iGallino et all (|1998l) and lArlandini etaLl (|1999T ) desig- 
nated the standard case for the mass of the 13 C pocket 
as that which best reproduced the S.S. main s-process 
component in low-mass (1.5 and 3.0 Mq) AGB models 
with [Fe/H] = —0.3. The s-process material in the S.S. 
was produced by a variety of AGB sources over many 
Gyr. Despite this fact, Figure [5] suggests that — at least 
for the elements with 56 < Z < 72 — the relative yields of 
low-metallicity, higher-mass AGB stars are not that dif- 
ferent from the more metal-rich, lower-mass AGB stars. 

In the M22 r + s group, 62% of the total amount 
of Z > 38 elements examined (excluding Ir and Th) 
originated in the s-process. The s-process contributes 
79% of the material to these same elements in the S.S. 
(jSneden et al.l 120081) . Hypothetically, if one wants to 
further enrich the heavy elements in the M22 r + s 
group to match the S.S. abundances, a greater fraction 
of s-process material (with respect to r-process material) 
needs to be added. In principle, then, these data support 
the general understanding that s-process enrichment oc- 
curs at later times than r-process enrichment. 



8. CONCLUSIONS 

One longstanding obstacle to properly interpret models 
of the s-process is having observations of pure s-process 
material outside the S.S. to compare with, especially 
since nearly all stars contain at least a trace of r-process 
material. Here we provide one solution to this problem 
by deriving the abundance patterns in two related groups 
of stars in the metal-poor GC M22. One group shows 
an r-process pattern with no detectable enrichment by 



s-process material (the r-only group), while the other 
group shows an additional s-process enhancement (the 
r + s group). By subtracting the r-process abundance 
pattern of the former from the r + s abundance pattern 
in the latter, we explicitly remove the r-process contri- 
bution to reveal the s-process "residual." 

The s-process abundance pattern in M22 strongly dis- 
favors low mass (M < 3 Mq), low-metallicity AGB mod- 
els. Although no published model results span the ap- 
propriate range of AGB masses at the metallicity of M22, 
the limited predictions available for more massive AGB 
stars at low metallicity fit the data better, especially the 
moderate Pb enhancement. Predictions for M = 4.5 and 
5.0 M© AGB models at [Fe/H] = -1.6 and -2.3 do fit the 
M22 s-process abundances, although 3 < M < 4.5 Mq 
models cannot be excluded because no predictions are 
available. The neutrons that fuel the s-process in these 
models mainly originate in the 22 Ne(a,n) 25 Mg reaction, 
which requires higher activation temperatures than the 
13 C(a,n) 16 reaction. In principle this could explain ob- 
served overabundances of K, Ca, Sc, V, Co, Ni, Cu, and 
Zn in the r + s group. We also calculate the r- and 
s-process fractions of each n-capture element. This ap- 
proach assumes nothing about the r- and s-process frac- 
tions in S.S. material. We encourage investigations of 
s-process nucleosynthesis in models with the appropriate 
metallicity and AGB mass range to better understand 
the origin of the heavy elements in M22. More generally, 
we hope that these data will serve as useful benchmarks 
for modeling and interpreting s-process abundance pat- 
terns and enrichment in low metallicity environments. 

Furthermore, these abundances can help interpret the 
enrichment history of M22. The Z > 27 abundance pat- 
tern in the M22 r-only and r + s groups bear striking re- 
semblance to the (physically unrelated) GCs M5 and M4, 
respectively. The r + s group in M22 may share a simi- 
lar enrichment history to M4 and possibly the metal-rich 
group in NGC 1851. If the s-process in M22 did originate 
in more massive AGB stars, this places strong constraints 
on the timescale for chemical enrichment, particularly in 
attempting to explain why the majority of metal-poor 
GCs do not show similar signatures of s-process enrich- 
ment. AGB models that can simultaneously explain the 
observed abundance patterns resulting from both proton- 
and neutron-capture reactions (the light clement varia- 
tions and s-process enrichment) should prove enlighten- 
ing in this regard. 
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APPENDIX 

LINE-BY-LINE MEAN OFFSETS 

We have calculated line-by-line mean offsets for n-capture species whose abundance is derived from three or more 
lines. Such information is useful when comparing abundances from different studies that use a small number of non- 
overlapping lines. In Table [9l we list the species (columns 1 and 6), wavelength (columns 2 and 7), average offset from 
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TABLE 9 
Line-by-line Mean Offsets 



Species 


A (A) 


(A) 




CJ 




Species 


A (A) 


(A) 




(J 




(Til 


Y II 


4883.68 


-0.009 


0. 


.081 


0.033 


Nd 


II 


4021. 


33 


-0.081 





.108 





.022 


Y II 


4982.13 


-0.003 


0. 


.086 


0.035 


Nd 


II 


4059. 


.95 


-0.095 





.144 





.029 


Y II 


5087.42 


0.005 


0. 


076 


0.031 


Nd 


II 


4232. 


37 


-0.034 





.074 





.015 


Y II 


5119.11 


0.007 


0. 


046 


0.019 


Nd 


II 


4446. 


38 


-0.013 





.051 





.010 


Y II 


5200.41 


-0.003 


0. 


.051 


0.021 


Nd 


II 


4462. 


98 


0.216 





.227 





.046 


Y II 


5205.73 


0.034 


0. 


045 


0.018 


Nd 


II 


4465. 


06 


0.054 





.067 





.014 


Y II 


5289.82 


-0.032 


0. 


.095 


0.039 


Nd 


II 


4465. 


59 


-0.006 





.051 





.010 


Zr II 


4050.33 


-0.200 


0. 


230 


0.163 


Nd 


II 


4501. 


81 


0.020 





.069 





.014 


Zr II 


4613.92 


0.145 


0. 


180 


0.127 


Nd 


II 


4567. 


61 


0.002 





.048 





.010 


Zr II 


5112.28 


0.055 


0. 


.066 


0.047 


Nd 


II 


4645. 


76 


-0.012 





.048 





.010 


La II 


3988.51 


-0.216 


0. 


.230 


0.064 


Nd 


II 


4706. 


.54 


0.068 





.075 





.015 


La II 


3995.74 


-0.115 


0. 


.185 


0.051 


Nd 


II 


4797. 


15 


-0.098 





.141 





.029 


La II 


4086.71 


-0.076 


0. 


.130 


0.036 


Nd 


II 


4825. 


.48 


0.007 





.041 





.008 


La II 


4322.50 


-0.085 


0. 


.098 


0.027 


Nd 


II 


4859. 


03 


0.023 





.042 





.009 


La II 


4662.50 


0.000 


0. 


031 


0.009 


Nd 


II 


4902. 


01 


0.084 





.091 





.019 


La II 


4748.73 


-0.083 


0. 


.093 


0.026 


Nd 


II 


4914. 


38 


-0.008 





.023 





.005 


La II 


4804.04 


0.068 


0. 


.074 


0.021 


Nd 


II 


5089. 


83 


=0.0 a 





.065 





.013 


La II 


4920.98 


0.161 


0. 


.172 


0.048 


Nd 


II 


5092. 


79 


-0.024 





.059 





.012 


La II 


4986.82 


0.034 


0. 


.071 


0.020 


Nd 


II 


5130. 


.59 


-0.059 





.073 





.015 


La II 


5114.56 


0.089 


0. 


.100 


0.028 


Nd 


II 


5132. 


33 


-0.055 





.102 





.021 


La II 


5290.84 


-0.074 


0. 


.097 


0.027 


Nd 


II 


5234. 


19 


-0.034 





.047 





.010 


La II 


5303.53 


0.068 


0. 


.092 


0.026 


Nd 


II 


5249. 


58 


0.070 





.083 





.017 


La II 


6262.29 


0.093 


0. 


.100 


0.028 


Nd 


II 


5255. 


51 


-0.006 





.038 





.008 


La II 


6390.48 


0.136 


0. 


.145 


0.040 


JMu 


II 


5293. 


16 


-0.052 





.090 





.018 


La II 


6774.27 


=0.0 a 


0. 


.075 


0.021 


In a 


II 


5319. 


81 


0.035 





.047 





.010 


Ce II 


4073.47 


-0.082 


0. 


.108 


0.028 


Sin 


II 


4318. 


93 


-0.034 





.067 





.024 


Ce II 


4083.22 


0.046 


0. 


.087 


0.023 


Sm 


II 


4434. 


.32 


0.053 





.083 





.029 


Ce II 


4120.83 


0.043 


0. 


.056 


0.014 


Sin 


II 


4467. 


.31 


-0.124 





.137 





.048 


Ce II 


4127.36 


-0.180 


0. 


.204 


0.053 


Sm 


II 


4536 


51 


0.038 





.055 





.019 


Ce II 


4137.65 


-0.053 


0. 


.146 


0.038 


Sm 


II 


4537. 


91 


-0.086 





.091 





.032 


Ce II 


4222.60 


0.027 


0. 


100 


0.026 


Sm 


II 


4591. 


81 


0.023 





.086 





.030 


Ce II 


4364.65 


-0.073 


0. 


100 


0.026 


Sm 


II 


4642. 


23 


0.075 





.089 





.031 


Ce II 


4418.78 


0.041 


0. 


056 


0.014 


Sm 


II 


4669. 


61 


-0.041 





.061 





.022 


Ce II 


4486.91 


-0.023 


0. 


.112 


0.029 


Sm 


11 


4719. 


84 


0.096 





.110 





.039 


Ce II 


4560.96 


0.016 


0. 


045 


0.012 


Eu 


II 


3907. 


11 


-0.246 





.252 


0. 


.178 


Ce II 


4562.36 


0.016 


0. 


046 


0.012 


Eu 


II 


4129. 


72 


0.009 





.087 





.062 


Ce II 


4572.28 


0.037 


0. 


075 


0.019 


Eu 


II 


6645. 


06 


0.237 





.258 





.182 


Ce II 


4582.50 


0.055 


0. 


095 


0.025 


Gd 


II 


4130. 


37 


0.173 





.190 





.134 


Ce II 


4628.16 


0.089 


0. 


102 


0.026 


Gd 


II 


4251. 


73 


-0.097 





.152 


0. 


.108 


Ce II 


5274.23 


0.032 


0. 


063 


0.016 


Gd 


II 


4498. 


29 


-0.076 





.127 





.090 


Ce II 


5330.56 


0.011 


0. 


087 


0.022 


Dy 


II 


3694. 


81 


-0.108 





.321 





.185 


Pr II 


4222.95 


0.012 


0. 


041 


0.024 


Dy 


II 


3983. 


65 


0.073 





.120 


0. 


.069 


Pr II 


4408.81 


-0.041 


0. 


050 


0.029 


Dy 


II 


4073. 


12 


-0.078 





.153 





.088 


Pr II 


5259.73 


0.012 


0. 


051 


0.030 


Dy 


II 


4449. 


70 


0.113 





.161 





.093 


Pr II 


5322.77 


0.017 


0. 


024 


0.014 





















Note. — For a given line, the mean offset (A) is computed as the average over all 
6 stars of the offset relative to the mean of all other lines of the same species in a given 
star. Exceptions are Gd II, whose mean is computed without 1-80, and Dy II, whose mean 
is computed without 1-27 and 1-80. 

a The log(g/) values for these two lines are not given in the literature and are derived here 
to empirically match the mean abundance derived from other lines in each of the 6 stars. 
See text for details. 

the mean abundance as derived for each of the 6 stars examined (columns 3 and 8), standard deviation of the average 
offset (columns 4 and 9), and standard deviation of the mean of the average offset (columns 5 and 10). 

Because the number of La n and Nd II lines examined is large, we have also derived empirical log(<?/) values for 
two lines not covered in the lLawler et a l. (2001a) and iDen Hartog et all (|2003l ) laboratory studies, La n 6774. 27 A 
(log(g/) = -1.77 ± 0.06) and Nd n 5089.83A (log(g/) = -1.27 ± 0.06). These lines are not used in determining 
the abundances in M22. This La II line i s often one of the only lines available in studies that target the red region 
of the spectrum. Uohnson fc Pilachowskil (|2010T ) provide empirical corrections to the abundance to account for the 
unknown HFS pattern of this line. In the M22 stars observed, the EWs of this line are all 10-20 mA, so the correction 
is approximately zero. 
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